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ABSTRACT 

A method is described for making estimates of the total 
emiasivity of hydrogen in the temperature and pressure ranges where 
hydrogen atoms predominate under equilibrium conditions. For a 
characteristic geometrical depth of 50 cm, and temperatures of the 
order of 12, 500°K and higher, with pressures of the order of 100 atmos 
and higher, the emissivity approaches unity (6 > 0. 95), while for 
temperatures of the order of 9500°K and lower, with pressures of the 
order of 10 atmos and lower, the emissivity approaches zero (fc^O.05). 
The variations of the emissivity between these approximate limits are 
shown graphically as functions of temperature and pressure with the 
geometrical depth set at 50 cm. The variation of the emissivity with 
geometrical depth is also shown graphically at 12, 600°K and 20 atmos. 
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I. 



INTRODUCTION 



A rocket utilising nuclear energy with hydrogen as working 
fluid has been considered in the past. It is therefore of come interest 
to obtain the emissivity of hydrogen atoms as a function of temperature 
and pressure in order to assess the importance of radiant heat trans- 
fer for high-temperature systems utilizing hydrogen. If the tempera- 
tures are sufficiently high, the hydrogen will exist almost completely 
in the atomic state, and hydrogen atoms will make the major 
contribution to the total emissivity. Furthermore, the treatment for 
hydrogen atoms provides a simple example of emissivity calculations 
for atomic systems. 

The emissivity, € , of a body is the ratio of the emitted radiant 
energy to the radiant energy emitted by a blackbody. The contributions 
to the emissivity of hydrogen atoms come from two sources: the con- 
tinuous spectrum and the discrete transitions. 

The blackbody curves may be obtained by plotting 
_ 2tt 

against oj for various temperatures. 

Here Rto dai (erg/cm^-sec) is the radiant energy emitted by a black- 
body, per unit area, per unit time, in the wave number interval 
into a solid angle of 2 tt steradians, T(*K) is the absolute temperature, 
C (cm/sec) is the velocity of light, ^(erg-sec) is Planck's constant, 
and (erg/°K) is the Boltzmann constant. If one considers the system 
of heated hydrogen atoms to be isothermal and uniform, then 
R^=R^ [_1- . This function can be plotted as a function of u) on 
the same graphs as Ru> • Here R^dco (erg/cm -sec) is the radiant 
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energy emitted by the hydrogen atoms, per unit area, per unit time, in 
the wave number interval doo, into a solid angle of 2 rr steradians, 

Pea (cm ^-atmos S is the spectral absorption coefficient, and 
X(cm-atmos) is the optical density. The spectral absorption coefficient 
Pto is set equal to Ro > where arises from discrete 

C 

transitions and from continuous spectrum contributions. To 

obtain the emissivity one takes the ratio of the actual area under the 

Koi vs. u) curve to the blackbody area, <r T , where 
-5 2 o 4 

<r* = 5.6699x10 erg/cm - K -sec. Examples of plots of Rto and 
P* as functions of 00 are reproduced in Figs. 1 to 5 for various 
temperatures at a representative pressure of 40 atmos. For conven- 
ience of presentation, R^/Rtu,„ox have been used as 

_ O O 

ordinates, where is the maximum value of Ru> . The con- 
siderations used for the calculation of Pc^ and hence will now 



be described. 



II. evaluation o?" the s-?fctp./i. arsor^tion coefficient, 



Pco , AND OF OPTICAL DENr^ITY, K. 
i-i. Calculation oi P<o 

The mass absorption coefficient in cm ' per gram of neutral 
hydrogen atoms for the continuous spectrum, X,/ • has been given by 
Unsold^ ^ ^ as 



1/'- fe4n^vrvc^° jy 






(1 - e"). 



(U 



Here vrv is the mass of the electron ( cjy is the charge of the 
electron ( e s o ), 111^ is the mass of the hydrogen atom 

= Jx c u) /'I? T , = R Kc/kl'* T , R is the Rydberg constant (cm’ ^ 

and yi is the principal quantum number of the hydrogen atom. To 

t C - 1 - 1 

convert Xv to , which has the units of cm -atmos , one must 

multiply by the atomic weight of hydrogen and divide by the absolute 
temperature and the molar gas constant, R , in cir. -atmos/rnolc- i . 

The contribution of the K*th spectral line whose center lies at 
cjo can generally be determined by computing the rris-sin.um value of 
* which is usually the value of R^(K)at the line center 
[- RI) (R)] multiplying this quantity by the intensity width defined 

as the wave number range for which R^j(K) R»,ox(N) • Detailed 
calculations have shov/n that this procedure will give only an error of 
approximately one percent for the contribution of the line to the total 
emissivity. i’or collected methods of computing the spectral absorp- 
tion coefficient of isolated lines see, for example, Penner and 
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( 2 ) 

Kavanagh . The maximum absorption coefficient for purely Doppler 

t 

broadened lines, P , is given by the equation 

P'- CS/cOoXm„ c/znXT)'" <2) 

where S (cm ^-atmos S is the integrated intensity of the spectral 
line under consideration. The integrated intensity, S , may be 
written as 

5- 5(MjN)[l-e>fpC-Aca)oAT]]^ 

3 

v/here Nl is the total number of atoms per cm per atmos, is the 

3 ^ 

number of atoms in the lower state per cm"^ per atmos. Here 5 is 
given by 

where f is the "oscillator strength", which has been tabulated by 
Unsold^ At equilibrium 

(--XccO,./.J^T)^ (5) 

where is the degeneracy of the lov/er state involved in the trans- 
ition, and Q is the partition function. In this case is taken 

as and = R(l- l/x) where n. is the principal quantum number 
of the lower state in the transition. In the present calculation the 
partition function is, for practical purposes, equal to the contribution 

from the ground state. For the temperatures and pressures considered 

(3) 

here it may be shown, using the estimates of Urey' ', that this 
procedure does not introduce errors greater than about one percent. It 
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is now possible to write Eq. (3) as 

S- 5'yu [i - cxpC-Acoj^/^T)] exp [Cl- 1 /f^X' (6) 

To estimate the intensity width, 2)0)4^- u7«| , at half the value 
of uses the well known curves of growth (see, for 

example, Unsold^ or Penner and Kavanagh^^^). Here is a 

wave number such that R^(b) - I^R„,a,CK) . The line shape para- 
meter is given by 

(7) 

where the natural half width has been neglected compared to the 
collision half width, . The Doppler half width, , is given by 

( 8 ) 

An estimate of the collision half width is obtained from the formula of 
Weisskopf^"^^: 

Jlr^~ c^yi^f/2.rrYncu)a, (9) 

3 

where >v is the number of molecules per cm . Combining Eqs. (4) 
and (9), one obtains 

4rc^S'Pr/lir^u>.^ ( 10 ) 

0 

where Pt is the pressure in atmos. Knowing cl and PX , it is 
possible to read from the curves of growth the value of the quantity ^ , 



(H) 
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where 

A = CJol R«. (12) 

Combining Sqs. (11) and (12), one obtains 

|u)4^,- cJo 1 = S’ Ao/ (13) 

If the Stark effect due to ionization is considered, may be 

taken as a first rough estimate, as Pj + PuT , where P(J is the 
absorption coefficient due to ionization broadening alone, and is 

the absorption coefficient due to broadening without considering 
ionization. A relation for the absorption coefficient due to ionization 
broadening is given by 

ot(AX)vb*ex pC'Jvcco./J^T) Wo (14) 

2 

where its Avogardo's number, and, <X ( a>\ ) (cm /atom) is the 

absorption coefficient per absorbing atom at a wave length aX ( A ) 
from l/tAo. The relation for (A X ) is 

0(CaX) = 4.2 1 C„ (a * (15) 

where N , the number of protons per cm , may be determined from 

( 5 ) 

the Saha equation, and has been tabulated by /\ller' ' and 

* For a more complete discussion of the Stark effect due to ionization 
see Ref. 5. Here broadening due to electrons is neglected. The 
electron effect has been discussed recently by B. Kivel, S. Bloom, 
and H. Margenau, Phya. Rev. 9^, 495 (1955). 
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Uneold^ K Using Eqs. (14) and (15), a check of ionisation broadening 
has been made at various temperatures and pressures using H« as 
an example. In all cases it was found that Ll‘ << I at 

the value of previously determined without considering Stark 

broadening, so that it seems justifiable to talce PuT = PuT except 
perhaps for the JLyman lines, where the factor exp j \J) is 

unity. 



The method described above to determine contributions to the 
emissivity due to discrete transitions does not apply in the case of the 
Lyman lines because they are too intense, and the blackbody curve 
along the intensity width of the lines varies appreciably in this area. 

At the temperatures considered here it is only necessary to determine 
modifications of the above method for Lyman , since the succeed- 
ing Lyman lines give negligible contributions because of the small value 
of the blackbody radiancy, Rto . at large bJ . When ( oJ - ) ie, 

for example, of ouch magnitude that it does not overlap an adjoining 
line, an expression for the discrete pressure broadening absorption 
coefficient, , has been given by Born^^^. For large values of 

\/(X, , Born's relation is 






(16) 



where 5= [(to- cJ»)/^g] (-^2.) , When ( CO “ b->o) iSi fo** example, of 
such magnitude that it overlaps an adjoining line, Born's method for 
pressure broadening and the treatment described above for the Stark 



^ - 



effect do not apply because one can no longer consider the contributions 
as resulting from "amall" perturbations of well defined energy states. 
The problem of the very broad Lyman line has caused difficulty in 

related astrophysical applications, and an exact treatment is not 
feasible at present. Therefore, a reasonable cut-off point has been 
chosen beyond which the above treatment io not justifiable. The cut-off 
point to the left of the line is taken arbitrarily as one half the distance 
between the Lyman ot and Lyman line centers. The cut-off point 
to the right of the line is unimportant since the blackbody contribution 
is negligibly small. The application of Lqs. (14), (15), and (16) to the 
Lyman lines results effectively in adding a sxrjall contribution of 
blackbody radiation from the left cut-off point to co » oo . It is felt 
that the neglect of the Lyman contribution from u) = o to the left 
cut-off point does not give appreciable error, but this approximation 
may well deserve further investigation. 

The contributions due to the higher series, Brackett, Pfund, 
etc. , have not been considered in the present computations. These 
series will give appreciable contributions only at much lower tempera- 
tures than were used. Cven if these contributions were blackbody from 
the left of the first Patchen line to OJ = O , the emissivity would in no 
case be increased by more than about 1 ^/o . 

B. Determination of Optical Density X 

The optical density, X , is defined as where is the 

partiad pressure of the hydrogen atoms (atmos) and Jl is the geometric 
depth (cm). The partial pressure is given as 
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(Ph= (17) 

where . % is the degree of dissociation of hydrogen molecules and Pr 
is the total pressure. The degree of dissociation, % , may be deter- 
mined from 

Kp = (4 ;c.yi - (IS) 

where Kp is the equilibrium constant for the dissociation of 
molecular hydrogen. Good estimates of Kp can be obtained by ex- 
trapolating data given in the National Bureau of Standards Tables^ 

Values of X- and are listed in Table I. From this table it is seen 

that hydrogen atoms may be expected to contribute significantly to the 
total emissivity. In a more exact treatment it would be necessary to 
consider the contributions of species other than hydrogen atoms which 
are present in relatively small amounts at the temperatures and 
pressures considered in the present investigations. 
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in. RE3U1-TS 

By use of the above procedures the total emissivity of hydrogen 
atoms in equilibrium with hydrogen molecules at various temperatures 
and pressures has been consputed. For most of the computations a 
characteristic geometrical depth of 50 cm has been used, and Figs. 1 
to 7 refer to that depth, as does the following discussion except where 
explicitly noted otherwise. In the ranges of ten.peratures and 
pressures chosen, the emissivity varies significantly with changes in 
the state variables. Typical plots used to determine the emissivity arc 
shown in Figs. 1 to 5, where the outer curve represents the blackbody 
curve, the cross-hatched area represents contributions from the 
continuous spectrum, and the solid area represents added contributions 
from the discrete transitions. The discontinuities in the continuous 
spectrum portions of Figs. 1 to 5 correspond to the onsets of the 
various continua; the Pfund continuum begins to contribute at cO = 4380, 
the Brackett continuum at uJ = 6850, the Paschen continuum at 
<x> = 12190, and the Balmer continuum at co = 27400. It is not 
necessary to consider the Byman continuum for the temperatures used 
here. As seen frons the figures, the continuous spectrum contributions 
to the emissivity increase as the temperature increases. The signifi- 
cant lines are located as follows: Paschen series, lJo = 5330, 7790, 
9130, 9940, and 10480; Balmer series, U)o = 15220, 20580, 23050, 
24400, 25150, and 25750; and Lyman series, - 82200. As may 

also be seen from the figures, the intensity width increases as the 
temiperature increases, and hence discrete transition contributions also 
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increase. In general, the continuous spectrum gives the major con- 
tribution to the cmissivity at high temperatures ( T > 920C°K), and the 
discrete transitions give the major contribution at low teinperatureo 
( T < 920C°t. ). In a similar manner it may be shov/n that for a constant 
temperature the emisci’/ity increases with pressure because the 
number density of emitters increases (fixed -H , X increasing) linearly 
with pressure; furthermore, the lines are broadened as the pressure 
is raised. 

The data listed in Table II are obtained from such plots as are 
shown in Figs. 1 to 5. These data are shown in graphical form in 
Figs. 6 and 7. From these figures it is seen that the cmissivity 
approaches unity (6 ^ 0.95) for temperatures of the order of 12, 500°K 
and higher, with pressures of the order of 100 atmos and Iiighor, while 
for temperatures of the order of 9500°K and lower, with pressures of 
the order of 10 atmos and lower, the cmissivity approaches zero 
(e < 0. 05). The graph in Fig. 8 shows the change in the total 
cmissivity as a function of geometrical depth for constant temperature 
(T= 12, 600*^K) and pressure ( Pr = 20 atmos). If either the tempera- 
ture or pressure is decreased, the cmissivity decreases and the 
€ vs ^ curve lies under the curve shown. Similarly, if either the 
temperature or pressure is increased, the 6 vs ^ curve lies above 



the curve shown. 
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TABLE II. EMI3SIVITY OF ATOMIC HYDROGEN AS A FUNCTION 
OF TEMPERATURE AND TOTAL PRESSURE 
( X = 50 cm). 





8400°K 


9200°K 


10080°K 


11300°K 


12600°K 


Pt (atmos) 


e 


6 


6 


e 


€ 


10 


.014 


. 037 


. C85 


. 26 


.48 


20 


.030 


.071 


. 16 


. 35 


.66 


40 


. 063 


. 14 


. 28 


.57 


.83 


70 


. 11 


. 24 


. 43 


.74 


.92 


100 


. 15 


. 32 


.55 


. 82 


.96 


150 


. 22 


.44 


.69 


.89 




200 


. 30 


.54 


.75 







I 

I 






-!'■)- 



£ 

3 




• ••••••••• 

^ O O O 0 ; 0;0 o o o 

• c 



10000 ZOOOO 30000 40000 50000 60000 70000 80000 90000 lOOOOO 

CO (cm') 

i{». 1. Plot for detemimng eimssivity at 8400 K and 40 atmos = 50 cm). 



t06 



u 




• • • • • • fl « 

^ OO O 0-0 5O o o o 

f c 



10000 20000 30000 40000 50000 60000 70000 80000 90000 100000 

W (CTT, ') 

Fig. 2 . Plot for determining emisBivity At 9»?00°K and 40 atmos = 50 cn ). 



90i 



H 

d 



r - 



o , 

3.2 

§ 

■2 S 

c 0 
o A 
O 



CQ 

n 

o 




10000 20000 30000 40000 50000 60000 70000 80000 90000 100000 

LO (cm"') 

Fig. 3. Plot for determining emissivity at 10,080^K and 40 atmos (X= 50 cm). 



i 

I 



0.90 



1 




10000 7.0000 30000 40000 50000 60000 70000 80000 90000 lOO'-Ou 

^ o n 

fig. 4. Plot for detern'.ining ei'r'iscivity at 11,300 K and 40 atmos 50 t_rr). 



t06 *0 



-19- 




10000 ?0000 30000 40000 .jOOOO 60000 70000 30000 90000 100000 

U3 (crn'') 

Fig. 3. Plot for determining ert issivity i t 1?, ^ 00°K ond 40 -tmos ( cn ). 



- 20 - 




. . . . . . . • . . vO 

,-hOOOO o ooo o 

ci 

MJ 



0 50 100 150 200 

{c.trrjos) 

Cmissivity of hydrogen aton s at v .r’ous ternpcratux es c.b a function of total pressure iJl,- 50 cm). 



1 



4 










I 




OSI = M t06 * 



- 21 - 




'4J 



P'lg. 7, ^missivity ot hydrogen atoms at various total pressures as a 1 uiction of temperature (j^= 50 cm). 



- 22 - 




20 40 f^O 80 100 120 140 160 180 200 

Z (cm) 

Fig. 8, Emissivity of hydrogen ritoms as a function of geometrical depth ( I = 12,600°K, P = 20 atmos). 



I 



♦ 



-23- 



APPENDIX 

A. Sample Calculation 

A detailed computation is presented here in order to show an 
application of the previously described procedure more clearly. This 
specific calculation is for the emissivity at lOOSO^K and Pj = 40 atmos 
with ^ = 50 cm. The Vcilues used for the various constants are as 
follows; 

c = 

M = 

a~ = 

yn = 

R = 

R' = 

1. Blackbody Curves , 

the relation = 2 Tr^c^ / ^xp For T = 10080*^K 

this becomes Rl,= 3.74x10 ^co ^[exp ( 1 . 427x1 0”^cO) - l3* Values of 

vs ui are shown in Table A-I. The area under the blackbody 
curve, <7- is 5.85x10^^ erg/cm^-feec. 

2. Contributions from the Continuous Spectrum. The mass 
absorption coefficient (Eq. (1)] at 10080®K is given by 



2.9979 X 10^^ cm/sec 
-27 

6.6238 X 10 erg-sec 
1. 3803 X 10"^^ erg/®K 
5.6699 X lO"^ erg/cm^-sec-°K 
4. 8022 X lO’^® esu 
9. 1072 X 10'^® g 
1.6819 X 10"^^ g 
1 . 0968 X 10® cm” ^ 

82. 06 cm®-atmos/mole-°K 

The blackbody curves are plotted using 




I 

I I 
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-%l = (?.6 2^ lo'^)[l- e%lp(rlA^^ - ^o'co)] G/coJ (A-1) 

where G~ 2 + € Vz-w-i . Values of G obtained from Unsold^^^ 

are as follows:* 





G 


4.38 to 6.85 


0. 0596 


6. 85 to 1 2. 19 


0. 101 


12.19 to 27.4 


0. 309 


27.4 to 109.7 


6. 43 


The spectral absorption coefficient, 


C 

Ru , is 



Pj= X'^ (/oo8)/R'T, (A-2) 

from which = 1.216xl0”^ . The optical density, X “ (f„ * 

is 1998 cm-atmos for this case ( " 39.95 atmos). The term 

1^1 - exp (- X )] can now be determined, and hence R^, can also 

C C 

be obtained. Table A-II lists values of Puj and Pco for various . 

The area of contributions from the continuous spectrum is found by 

11 2 

graphical methods to be 1. 118x10 erg/cm -sec. 

3. Contributions from the Discrete Transitions. To calculate 
p' as given in Eq. (2), it is first necessary to compute 5 [Eq. (4)] 



♦ Unsold does not list values of G for 9200*^K and 11300°K; however, 
these can easily be calculated as indicated by the formula for G . 
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and 5 [Eq. (6)], The value of N (atoms /crn^-atmos) in Eq. (4) 
can be obtained from the perfect gas law, PtV = T, where n is 
the number of moles and v ie the volume (cm^). Rewriting this 
with R = No , one obtains N= - j/J?T . Hence Sq. (4) 

becomes 

5'= 6.45 X 

where f is taken from the tabulated values of Unsold^ The 
relations for S for the Balmer and Paschen lines, respectively, 
become from Eq. (6), 

5 a.ofZ- exp(-f,42 7 X ( 6 * 0 J^] s', 

and (A-4) 

5 = 8.046'^ - exp(-|.427x IO‘"^tUo)] 5. 

Equation (2) is then written 

p' = /.3/5 X 10“^ S/cO^. (a-5) 



At the given temperature and pressure 
.4 < , 2 



= 3.58x10'^ tdo [Eq. (8)], 
and <X = 4.71x10^ 5 / oof [Eqs. (7), (8), and (9)J . The terms CU 
and logjQ 10.6 PX can now be determined and ^ then found 
from the curves of growth. The intensity v/idth from Eq. (13) is now 



2|o)j^-U)o| = 2.40 3 (A-6) 

Table A-III gives the necessary data to obtain 2 for the 

spectral lines considered, and the last column of this table gives the 
value of 2 computed from these data. The additional area 
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contributed by each line is usually obtained by multiplying dJoj 

by (Rft> " Ruj ) at 6 l>o . If 2|c0|^- U)<,| is so large that either 

_ C 

or deviate appreciably from linear variation along the entire 

intensity width, then the area must be found by graphical methods. These 
additional areas are listed below: 



Line 

3 

4 

5 

6 

7 

8 



Balmer 

2 

Area (erg/cm -sec) 
3. 13 X 10^^ 

0. 650 
0.224 
0. 098 
0. 056 
0.037 



Line 

4 

5 

6 

7 

8 



Paschen 

2 

Area (erg/cm -sec) 
0. 635 X 10^° 

0. 204 
0. 125 
0. 052 
0. 037 



The Lyman contribution is taken as blackbody from the left cut-off 

4 

point (CO =7. 46x10 ) to cO = oo . This area is found to be 
10 2 

0.250x10 erg/cm -sec. 

Hence, the total area contributing to the emissivity, i. e. , the 
sum of the continuous spectrum and discrete transition contributions, 
is 16.43x10^^ erg/cm^-sec , and therefore € = 1.643x10^^/5.85x10^^ = 



0. 281. 
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TABLE A-I. R.1 AS A FUNCTION OF CO ( T = lOOSo'^K). 



<^/ 10 ^ 


R«/10^ 


1 


0 . 244 


3 


1.90 


4 . 38 


3.63 


6 . 85 


7.30 • 


10 


11.80 


12 . 19 


14 . 40 


15 


16.83 


20 


18.22 


25 


16.94 


27 . 4 


15.70 


30 


14.02 


32 . 5 


12.44 


35 


11.05 


40 


7.87 


45 


5.50 


50 


3.74 


60 


1.55 


70 


0.589 


80 


0 . 212 


90 


0.072 



100 



0 . 025 
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TABLE A-II. 


II 


AND AS FUNCTION OF 

40 ATMOS, X = 50 CM). 


UJ ( T = 10080°K, 


60/10^ 




Pj/10"'^ 




00 

• 




3.82 


1. 94 


6.85 




1.35 


1.73 


6.85 




2. 29 


2. 68 


10 




0.900 


1.94 


12. 19 




0.540 


1.47 


12. 19 




1.65 


4.05 


15 




0.944 


2.90 


20 




0.426 


1.49 


25 




0.225 


0.748 


27.4 




0. 172 


0.531 


27. 4 




3.58 


8.05 


30 




2.75 


5.94 


32. 5 




2. 18 


4. 40 


35 




1.75 


3. 25 


40 




1.17 


1.65 


45 




0.826 


0.839 


50 




0.600 


0. 423 


60 




0. 348 


0. 104 


70 




0. 220 


0. 025 



TABLE A-m. INTENSITY WIDTH AND MISCELLANEOUS DATA FOR COMPUTING INTENSITY 
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